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Abstract 

We consider the radiative transition — > / 7, which is a sensitive probe 
of the nature of the / (980) particle. Using the QCD sum-rule technique, 
we estimate the branching ratio of such decay mode to be: B((f> — > /07) = 
(2.7 ± 1.1) 10~ 4 , in fair agreement with present experimental data. As for 
the structure of the fo, the result suggests a sizeable ss component; however, 
this result does not exclude the possibility of further components and allows 
a more complex structure than indicated by the naive quark model. 



The quark model provides a rather good description of hadrons, which fit into suitable 
multiplets reasonably well. In its simplest version the model then interprets mesons as 
pure qq states. Scalar mesons present a remarkable exception to this successful scheme. 
Indeed, the nature of these mesons is not established yet Q. There are more scalars 
than can fit into one quark model multiplet. Consequently, some of these states could 
be either glueballs or admixtures of quark and gluonic states, or belong to multiquark 
multiplets. A particular feature of some of these particles is that they appear to be rather 
wide 0, |, f|, They have very short lifetimes and large couplings to hadronic channels, 
such as KK or titt. This might suggest that they can be identified as composite systems 
of hadrons, or that they spend an appreciable part of their lifetimes as such states. This 
could be the result of hadronic dressing, whereby the strong interaction enriches a qq state 
with other components such as KK^, \ttt]), etc. Such a viewpoint could also explain why 
the scalar mesons seem to contradict the OZI rule. Since the two mesons composing the 
state in which they spend much of their lifetime may readily annihilate to qq, leading to 
a subsequent OZI allowed decay. 

In this letter, we focus on the structure of the /o(980) and the possibility of gleaning 
information about this from radiative <p decays. According to the quark model, the /o(980) 
should be an ss state, an interpretation supported in Refs. || [7|, |J. However, this does 
not explain its mass degeneracy with the ao(980), that should be a (uu — dd)/y2 state. 
There are also suggestions that the /o(980) could be a four quark qqqq state ||. In this 
case, it could either be nucleon-like fllP] , i-e. a bound state of quarks with symbolic 
quark structure ss(uu + dd)/y/2^, or deuteron-like, i.e. a bound state of hadrons, which 
is usually referred to as a KK molecule [§, [II], |12|, In the former of these two 

possibilities, the mesons are treated as point-like objects, while in the latter they should 
be considered as extended objects. Some objections have been raised against the KK 
molecular model P, Ibl. In particular, such an interpretation requires a width smaller 



than the binding energy of the molecule itself, which has been estimated to be e ~ 10 — 20 
MeV [Q], in contrast to the measured width lying in the range 40 — 100 MeV |]15|] . 

Various ways have been suggested of clarifying the situation, such as the analysis of 

T(<f> -> a 7) 



the /o ~~ > 77 decay [16, O] or of the ratio — — - Jill. In the naive quark model, 

r(0 -> /07) 

for example, it is expected that B(<f> — > /07) and B(<p — > a 7) would differ by a factor of 



10. Moreover the rate for — > /o7 m ay distinguish among the different possibilities [11 



1 Within the same framework the isovector partner ao(980) is written as ss(uu + dd)/V2. 

2 



since, according to the existing theoretical estimates, the expected branching ratio would 
be as high as 10~ 4 in the qqqq case, (9(10~ 5 ) in the ss case. For a KK molecule, the 
branching ratio clearly depends on its size. For a compact state this is ~ 7 • 10 -5 , while 
for a diffuse, deuteron-like system, it is down below 10~ 5 



ff 



From the experimental point of view, the PDG value []15 



£(0^/o7) = (3.4 ±0.4) 10" 



(1) 



stems from averaging the results of the CMD2 [|T^j and SND [19| collaborations, analysing 
7T + 7r _ 7, 7r°7r°7 and 57 final states. What is more, a significant improvement is expected 



at the factory DA$NE [20], where the first results give: 

#0 -> /o7 -> 7r°7r°7) = (0.81 ± 0.09(stat) ± 0.06(syst)) x 10~ 4 



and 



at 90% C.L. [BTJ 



B(<P -> /o7 -> vr"7T + 7) < 1.64 x IO- 



CS) 



(3) 



The present letter is devoted to analysis of the radiative decay — > /07 using QCD 
sum-rules [^2|, which we previously applied to the radiative transitions to 77,77' p3| . 
That the /o(980) couples significantly through ss components has long been known [] 



from its appearance as a peak in J /if) — > <f>f [^5| and D s — > 7r/ [p£|, as discussed in 
Refs. 0, and in more detail in [^] . Our calculation relies on the assumed coupling of 
the /o to the scalar ss density. As a preliminary, we evaluate the strength of this coupling 
using two point QCD sum-rules. The result will then be exploited in the three point QCD 
sum-rule evaluation of the relevant quantity needed to compute B(<p — ► /07). 



The coupling of the /o(980) to the scalar current J s = 
of a constant /: 

(0| J s |/ (p)> = m f0 f . 



ss can be parametrized in terms 



(4) 



In order to compute this parameter by QCD sum-rules, we consider the two-point corre- 
lator: 



T{q 2 



d A xe iq - x (0|T[J s (x)J st (0)] |0) 
2 And noticed more recently by Delbourgo et al. (24j for (f> — > /07. 



(5) 
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which is given by the dispersive representation: 

\ roc p(s) 

T(q 2 ) = — ds + subtractions . (6) 

7T J4m1 s — q 2 

In the region of low values of s, the physical spectral density contains a 5— function term 
corresponding, in the small width approximation, to the coupling of the /o to the scalar 
current. Picking up this contribution and dropping possible subtractions which we discuss 
later, we can write: 

T{q > ) = ^JL + Lr ds tM, (7) 

m f — Q 71 J s o s — q 2 
assuming that the contribution of higher resonances and continuum of states start from 
an effective threshold sq. On the other hand, the correlator T(q 2 ) can be computed in 
QCD for large Euclidean values of q 2 , by using the Operator Product Expansion (OPE) 
to expand the T-product in Eq. (f|) as the sum of a perturbative contribution plus non- 
perturbative terms which are proportional to vacuum expectation values of quark and 
gluon gauge-invariant operators of increasing dimension, the so called vacuum conden- 
sates. In practice, only a few condensates are included, the most important contributions 
coming from the dimension 3 < qq > and dimension 5 < qgaGq >. 

In the QCD expression for the two-point correlator considered, the perturbative term 
can also be written dispersively, so that: 

oert^ 

T 

where the spectral function (f ert and the coefficients d 3 , d§ can be computed in QCD. The 
next step consists in assuming quark-hadron duality, which amounts to the claim that the 
physical and the perturbative spectral densities give the same result when integrated 
appropriately above some Sq. This leads to the sum-rule: 



1 roo rp ert (^\ 
T QCD {q 2 ) = - ds P - { Ar + d 3 <ss> +d 5 < sgaGs > +... 

7T J Ami S — q 2 



m% P ^ 1 p ^ (f er \s) 
m % ~ Q 2 71 Jim i s — q 



rso ff ert (s) 

/ ds -J- + d 3 <ss > +d 5 < sgaGs > + . . . (9) 

J Ami S — q 2 



This expression can be improved by applying to both sides of Eq. (|9]) a Borel transform, 
defined as follows: 

B\HQ 2 )] = Um n , _ ^ M2 j^i-QT " HQ 2 ) , (io) 



n- I)!' ^ ' \dQ 2 



4 



where T is a generic function of Q 2 = —q 2 . The application of such a procedure to the 
sum-rules amounts to exploiting the following result: 



B 



1 



[s + Q 



2\n 



exp(-s/M 2 



(M 2 ) n (n- 1)! 

where M 2 is known as the Borel parameter. This operation improves the convergence of 
the series in the OPE by factorials in n and, for suitably chosen values of M 2 , enhances 
the contribution of low lying states. Moreover, since the Borel transform of a polynomial 
vanishes, it is correct to neglect subtraction terms in Eq. (|^), which are polynomials in 
q 2 . The final sum-rule reads: 



m 



fo 



rexp 



m 



M 2 



+ m s exp I - — 2 



<ss> 



( n m 2 mt \ ^ 1 / 



mz 



2M 2 I 



(12) 



In the numerical evaluation of Eq. fll2|) we use < ss >= 0.8 < qq >, < qq >= (— 0.24) 3 
GeV 3 , < sgaGs >= 0.8 GeV 2 < ss >, m/ = 0.980 GeV. The strange quark mass 
is chosen in the range m s = 0.125 — 0.160 GeV, obtained in the same QCD sum-rule 
framework |28]. The threshold is chosen below a possible /o(1370) pole and varied between 
So = 1.6 — 1.7 GeV 2 . Since the Borel parameter has no physical meaning, we look for a 
range of its values ("stability window") where the sum-rule is almost independent on M 2 . 
Such a window is usually sought in a restricted interval of values of the Borel parameter 
chosen by requiring that the perturbative contribution is at least 20% of the continuum 
and additionally requiring that the perturbative term is greater than the non-perturbative 
contribution. The stability window for M 2 is selected in [1.2,2] GeV 2 , as seen in Fig. [I], 
where, taking into account the uncertainty on m s , we obtain the coupling: 



/= (0.180 ±0.015) GeV . 

This result will be used in the analysis of the decay <fi - 
following. 

The relevant matrix element describing the transition 
vector current J M = s^^s, can be parameterized as follows: 



(13) 

/o7 as we shall see in the 
— > f induced by a strange 



l/o(g 2 )| J M ei)> = F 1 (q 2 )(q 1 -q 2 )e^ + F 2 (q 2 )(e 1 -q 2 )(q 1 + q 2 ) 

+ Fsiq 2 ) (ei-q 2 ) q» 



(14) 
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Figure 1: Coupling of the /o to the scalar current as a function of the Borel parameter M, 
for m s = 0.140 GeV. The solid curve corresponds to the higher threshold Sq = 1.7 GeV 2 , 
the dashed curve corresponds to s = 1.6 GeV 2 . 



where q = q% — q%. In order to consider the radiative decay 
amplitude 

1 



.A(0(gi,ei) -> fo(q2h(q,e)) 



— > /07, one needs the 
e e*^ [Fi(0) (q r q 2 ) e lM + F 2 (0) (e r g 2 ) (gi+g 2 )J , (15) 



where the charge of the strange quark has been explicitly written. Eq. (JTSI) shows that 
only two of the three form factors appearing in Eq. flT4]) are actually needed. Furthermore, 



gauge invariance requires that q^ • [-^1(0) (q% ■ q^)^ + -^(0) (ei • q%) (li + — 0, which 
relates the values of Fi and F 2 at q 2 = 0: 



F 2 (0) = ^(0) 



2(mJ - m\) ' 

In terms of -P\(0), the rate for the process we consider becomes: 



(16) 



r(</> - /07) = « [*i(o)] 



2 (m 2 -m) Q ){ml + m 2 Q ) 2 



216 mi; 



(17) 



Three-point QCD sum-rules can be applied to evaluate the form factor Fi(q 2 ). We consider 
the three-point function: 



rWg 2 , q^ q 2 ) = i 2 / (fx d 4 y e~^ x e^ y (0| T[J%) J„(0) J„(x)] |0) 



where J s has been defined above and J u = 37^5 is the vector current. The correlator 
Eq. (|L8|) can be written in terms of invariant structures as follows: 

IV(gi, qj, q 2 ) = U(ql, q%, q 2 ) + n^g 2 , q 2 , q 2 ) q Xil q Xv + ■■■ (19) 

and a QCD sum-rule can be built up for the structure II(g 2 ,g|,g 2 ). The method closely 
follows the one described for the two-point sum-rule. We assume IT(g 2 ,g|,g 2 ) obeys a 
dispersion relation in both the variables q\,q\: 

n(fl?, ql f) = -J dsi / ds 2 P^lft , (20) 
vr 2 J J (s x - q{){s 2 - qi) 

with possible subtractions. Such a representation is true at each order in perturbation 
theory and, as is standard in QCD sum rule analyses, it is assumed to hold in general. 
In this case the spectral function contains, for low values of s x , s 2 , a double 5— function 
corresponding to the transition — > / - Extracting this contribution, we can write: 



n / n 2 2 jx ™/o / m U Fiiq 2 )^ ■ q 2 ) 1 /- / ad (s l5 s 2 , g 2 ) 

Ail9n92>9J- 7^2 „2W_2 ~\ V ~\, " S l / " S 2 

^7o 



(raj - - <&) r *" 2 h "° 2 (si - q 2 )(s 2 - qi) ' 

(21) 

where subtractions are neglected as later they will vanish on taking a Borel transform. 
The parameter / appearing in the previous equation is just the coupling of the f to the 
scalar current, computed previously. Deriving an OPE-based QCD expansion for II for 
large and negative qf, ql and q 2 , one can write: 

oert I „ „ „2\ 



rr/ 2 2 2\ 1 A A QWJ l / - ^ , ^ ^ 

1%!, q 2 , q ) = n / asi / «s 2 - 2T7 2T + c 3 < ss > +c 5 < sgcrGs > +... . 

Tl 2 JAmj J Ami {s X - q()(s 2 - qi) 

(22) 

Invoking quark-hadron global duality as before, we arrive at the sum-rule: 

m f0 f m U F x {q 2 ) {q 2 - m 2 - m 2 J _ 1_ f ^ ^ p^\ Sl , s 2 , q 2 ) 



2 ( m l - Qi)(™% - ql) k 2 Jd (si - qf)(s 2 - ql) 

+ C3 < ss > +c 5 < sgaGs > +... , (23) 

where the domain D should now also satisfy the kinematical constraints specified below. 
After a double Borel transform in the variables — q\ and — q 2 , we obtain: 



m f Jm^f^(q 2 - mj - mjj F x (g 2 ) exp ( - 



^L dsids2exv (-w 2 



M 2 Mf 



- — - — ^( Sl)S2 ) (24) 
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+ exp 



m: 



M 2 



m s 

Ml 



<ss> 

m 2 s q 2 (2m 2 s 



2 2 

g 2 + 2m 2 - s ; 



+ 



+ < sgaGs > 



mz 



1 q 
3 + 3M X 2 



+ 



2M{M 2 2 
2q 2 + m; 



M{ 



3M? 



g 2 (5m 2 — 2q 



where: 



3m, 



{(4m 2 + si-s 2 + g 2 )[( 



6M/M 2 2 

si + s 2 - g 2 

si + s 2 



2 2 
AMn 



J_ _1_ 

Mf + M 2 4 



m 2 (3g 2 + m\ 



12Mf 



4sis 2 + 4g sis 2 

3/2 



,2\2 



4sis 2 



(25) 



The integration domain D over the variables si, s 2 depends on the value of g 2 . For (— g 2 ) > 
s — 4m 2 , D is specified by: (s 2 )_ < s 2 < s 4m 2 < s x < s ; while, for (— g 2 ) < s — 4m 2 
D is bounded by: (s 2 )_ < s 2 < (s 2 )+ if 4m 2 < s\ < and (s 2 )_ < s 2 < sq if < 



s i < s , with: (s 2 )± = 2m 2 g 2 + (2m 2 — g 2 )si ± y S\q 2 (q 2 — 4m 2 )(si — 4m 2 ) /2m 2 and 
(si)-t- = 2m 2 g 2 + (2m 2 - g 2 )s ± y / s g 2 (g 2 - Am 2 s )(s - Am 2 ) /2m 2 . 

Since we consider the form-factor Fi(q 2 ) for arbitrary negative values of g 2 , we could 
perform a double Borel transform in the two variables Q 2 = —q\ and Q\ = — gf , which 
allows us to remove single poles in the Si and s 2 channels from the sum-rule. Our pro- 
cedure is therefore to compute the form-factor F\(q 2 ) and then to extrapolate the result 



to g 2 = 0. In the numerical analysis we use: m^ = 1.02 GeV, 



0.234 GeV (obtained 



from the experimental datum on the decay to e + e 



We compute the result for two 



values of the threshold: s f = 1.8, 1.9 GeV 2 . s coincides with the f threshold chosen 
as for the two point function. The extrapolation to g 2 = shown in Fig. |] gives: 

F l (0) = 0.34 ±0.07 , (26) 



which, using T(0) = 4.458 MeV flHJ and Eq. (|17D, gives: 

B(0-/„7) = (2.7±l.l) 10 



(27) 

Both the results Eq. flllf ) and Eq. (|2~6|) have been derived without the inclusion of radiative 
a s corrections, an approximation which is usually believed more accurate for the three 
point sum rule, where the 0(a s ) corrections are expected to cancel in the ratio of a 
three-point and a two-point function. 
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Figure 2: Form factor Fi(q 2 ). The dashed and dotted lines are the highest and the lowest 
curves obtained varying the set of parameters entering in the sum rule (p4|). The isolated 
point on the right is the result of an extrapolation. The result Q2"E| ) corresponds to the 
central point on the right obtained by extrapolating the solid curve. 



Our result of Eq. ( p7|) is in reasonable agreement with the outcome of refs. [^, |TT 



where the decay is supposed to proceed through the chain <fi — > KK^f — > /07, and so 
depends on the coupling gf oK x 0- Their results are: B(<f> — > /07) = 1.9 x 1CT 4 [^9| and 
B((f> — > /07) = 1.35 x 1CT 4 PH . On the other hand, QCD spectral sum rules are exploited 
in ref. Q to predict B(<j> -> /07) = 1.3 x 10" 4 . 



A different strategy is proposed in [16], where the experimental datum is assumed 



together with the structure /o(980) = nncosO + sssin#, where nn = (uu + dct)/y2 and 
9 is a mixing angle. A theoretical prediction is derived describing the particles (<f), /o) 
through wave functions depending on the radii of the mesons. Such a prediction is then 
compared to the experimental datum in order to constrain the mixing angle. 

Although our result of Eq. fl2"T|) is affected by a rather large uncertainty, it is in 



agreement with the available data [|T^, [19| . Since our sum rule analysis is based on the 
hypothesis that the / (980) couples to the scalar ss current, this agreement leads to the 
conclusion that an ss component is present in such a state. However, our branching 
ratio is an order of magnitude larger than the naive quark model gives for a pure ss 



3 Such a coupling is taken from ref. 



state. Our result is consequently consistent with the view that the / (980) is a meson 
with a basic qq composition, which spends a sizeable part of its lifetime in a two meson 
state, such as KK. This is in keeping with the analyses of |6|, |8], |32| that attribute 
such multi-hadron components to dressing. While the effect of KK couplings have been 
studied phenomenologically in Ref. |16| in a range of hadronic reactions, they have been 
dynamically calculated by Marco et al. [33[ explicitly for the radiative decay we study 
here and found to give £>(</> — > f j) = 2.4 x 10~ 4 , in reassuringly good agreement with 
our sum- rule result, Eq. (p7|). 
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